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Abstract. For a function V : R — > M that is integrable and compactly supported, 
we prove the norm resolvent convergence, as e — > 0, of a family S E of one-dimensional 
Schrodinger operators on the line of the form 



then the functions e~ 2 V(x/e) converge in the sense of distributions as e — > to 5'(x), 
and the limit Sq of S £ might be considered as a 'physically motivated' interpretation of 
the one-dimensional Schrodinger operator with potential S' . In 1985, Seba claimed that 
the limit operator So is the direct sum of the free Schrodinger operators on positive 
and negative semi-axes subject to the Dirichlet condition at x = 0, which suggested 
that in dimension 1 there is no non-trivial Hamiltonian with potential S'. In this 
paper, we show that in fact So essentially depends on V: although the above results 
are true generically, in the exceptional (or 'resonant') case, the limit So is non-trivial 
and is determined by the properties of an auxiliary Sturm-Liouville spectral problem 
associated with V. We then set U(£) = a^^) with a fixed \E' and show that there 
exists a countable set of resonances {ak}f?L_ 00 for which a partial transmission of the 
wave package occurs for So- 



1. Introduction 

Assume that V is an integrable function of compact support contained in the 
interval [—1, 1] and consider Schrodinger operators S £ on the real line given by 



Here e is a positive parameter, and one of the questions of our primary interest in this 
paper is the behaviour of S £ as e tends to zero. 




If the potential V satisfies the conditions 
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The motivation for this question stems from the fact that if V has zero mean and 
its first moment L £V(£) d£ is —1, then the functions e~ 2 V(x/e) converge in the sense 
of distributions as e — > to the derivative 5' of the Dirac delta-function. Therefore, 
if the Hamiltonians S e converge (in some topology) as e — > to a limit So, then it is 
natural to regard So as a realization of the Schrodinger operator with a potential 5'. 

Schrddinger operators in M. d , d > 1, with singular distributional potentials 
supported on a discrete set (such potentials are usually termed "point interactions") 
have attracted considerable attention both in the physical and mathematical literature 
over several past decades. One of the reasons for this is that such singular Hamiltonians 
have widely been used in quantum mechanics to model interactions in particles and 
atoms. For instance, as back as in 1931 Kronig and Penney [21] already suggested their 
model of electrons moving in the crystal lattices that used one- dimensional Schrodinger 
operators with periodic potentials of the form ot^2 neZ 8(- — n) as the corresponding 
Hamiltonians. Another reason is that Schrodinger operators with point interactions 
often form "solvable" models in the sense that the resolvents and other objects related to 
such operators can explicitly be calculated (see the books by Albeverio, Gesztesy, H0egh- 
Krohn, and Holden [2] and by Albeverio and Kurasov [5] discussing point interactions 
and more general singular perturbations of the free Schrodinger operators in M d and the 
extensive bibliography lists therein). This in turn allows one to discover new, unusual 
effects not seen for regular potentials as was the case, e.g. with the singular Wannier- 
Stark systems |14j . 

One of the most natural ways to define a Hamiltonian corresponding to a point 
interaction supported by x = is to first restrict the free Schrodinger operator onto the 
set of smooth functions that vanish in a neighbourhood of x = and then take a self- 
adjoint extension of the resulting symmetric operator S' . This approach was suggested 
for the first time by Berezin and Faddeev [8] in 1961 for three-dimensional Schrodinger 
operators with a potential a5; the operator S then has deficiency indices (1,1), and 
all self-adjoint extensions of S' Q form a one-parameter family i? 7 . In [8], the authors 
suggested a physically motivated choice of 7 for a given a based on the renormalization 
technique pertinent to the physical setting of the problem. The existence of nontrivial 
limits of Hamiltonians as the singular potential ad~ is approximated by a regular sequence 
of compactly supported ones was also demonstrated in [18]. We also notice that 
differential operators with singular coefficients can be studied within the framework 
of the contemporary theory of new generalized functions such as the Colombeau or 
Egorov's algebras, in which multiplication is well defined and which contain rich sets of 
(different) (^-functions obtained via approximating sequences with different profiles. This 
approach was applied in [7] to study the Schrodinger operator in M 3 with 5-potentials. 

In dimension 1, the symmetric operator S' has deficiency indices (2,2); thus the 
set of all its self-adjoint extensions forms a four-parameter family, and the problem of 
choosing a single Hamiltonian corresponding to a particular singular point interaction 
becomes more subtle. There are two families of extensions that have been studied 
especially well and have a clear physical interpretation. 
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The first one, denoted by S a j, is given by S a> sf = —f" on the domain 
dom^ := {/ G W*(R \ {0}) | /(0+) = /(0-) =: /(0), /' (0+) - /'(0-) = a/(0)} 
and corresponds to the Hamiltonian 
d2 

That such an identification is most natural is seen from the fact that various known 
approaches to definition of the above Hamiltonian (e.g. via the form sum pEJCIH], 
generalized sum 0, [10], approximation by regular potentials [9], and regularization 
by quasi-derivatives [28]) lead to the same operator S a ^s- In particular, if V n are 
regular functions converging weakly to the measure a8, then the corresponding regular 
Schrddinger operators 

converge in the norm resolvent sense to S a> s [S], [2J Ch. 1.3.2]. 

The second well-studied family, denoted Sp t $r, °f self-adjoint extensions of the 
operator S' in one dimension is given by Sp^rf = — /" on the set of functions 

domSp, s> = {/ G W?(R\{0}) | /' (0+) = f(O-) =: f (0), /(0+) - /(0-) = (3f'(0)} 

and is widely accepted as a model for Schrddinger operators with 8'- interactions (not 
(^'-potentials!), i.e. for the operators 

~ + P(-,6')S'. 

Here (/, ip) denotes the action ip{f) of a distribution ip on a test function /. Physically, 
this model corresponds to an idealized dipole of zero range at x — 0. Similarly to the 
case of 5-potentials, the operator Sp^i can be given as the limit in the strong resolvent 
sense of regular rank-1 perturbations of the free Schrddinger operator by f3(-,V n )L 2 V n 
whenever V n converge to 8' in the distributional sense [2JE]- Recently it was realized 
that Sp t s' is the limit in the norm resolvent sense of the operators under the special 
choice of K 0H5]. 

We note that one can also define Schrddinger operators with (^'-interactions 
supported by subsets T of R. Spectral properties of the corresponding Hamiltonians 
for discrete T were studied e.g. in [TTJED] and for T a Cantor set in [21]. Also, there are 
papers where more general rank-1 perturbations of abstract positive Hamiltonians H 
are considered, namely perturbations by a(-,ip)il), where ip is a generalized function in 
the domain of the operator (H Q + I)~ n , for n G N. The efficient way to treat such 
singular perturbations is via Krein's resolvent identity, see [3] and especially the papers 
by Kurasov [22] and Nizhnik [25] and the references therein. 

The situation is more obscure with definition of the Schrddinger operator with a 
potential 8', which might be considered as an alternative way of modelling idealized 
dipole of zero range at x = 0. Mathematically, such a potential should be interpreted 
as an operator of multiplication by the distribution 8'. The product S'f for smooth 
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enough / is a distribution, so that any Schrodinger operator chosen to serve as a putative 
model for such a physical system should be understood in the distributional sense. A 
natural approach to defining the operator is to approximate 8' by regular potentials (e.g. 
in the distributional sense) and then to investigate the convergence of the corresponding 
family of regular Schrodinger operators. As we explained above, one possible way to do 
this is to consider the family S £ of ([I]) with V of zero mean and nonzero first moment. 

This approach was realized by Seba [27]. In fact, he considered a family of 
Schrodinger operators with more general short-range potentials e~ a V (x / e) in theorem 4 
of [27] and claimed in particular that for a > | the operators S £ converge in the 
norm resolvent sense as e — > to the direct sum S- © S + of the unperturbed half- 
line Schrodinger operators subject to the Dirichlet boundary conditions at x — 0. 
One therefore could suggest that, in dimension one, no non-trivial interpretation of 8' 
potentials is possible. However, a careful analysis of the arguments given in [27] reveals 
that, although the proof is solid when a < 2, the term e 2 ~ a T e (k) in the expansion of 
the resolvent (S £ — k 2 )~ l does not vanish in the limit as e — > if a > 2, thus putting 
under question the validity of the result in this case. 

This was one of the motivations for us to re-examine the above convergence result. 
In this paper, we confirm part of the statement from [27] that for an arbitrary real-valued 
integrable potential V of compact support the family (JTJ of Schrodinger operators S £ 
converges in the norm resolvent sense to a limit denoted Sq. However, the limit So 
turns out to heavily depend on weather or not V is resonant. To define this notion, we 
consider the Sturm-Liouville operator M given by 

My = -y" + Vy (3) 

on the interval [—1, 1] subject to the Neumann boundary conditions y'(—l) = y'(l) = 
and call the potential V resonant if M has a non-trivial null-space and non-resonant 
otherwise. We show that in the non-resonant case, the limiting operator So is just the 
direct sum S_ ®S + of the free Schrodinger operators on R_ and R + respectively subject 
to the Dirichlet boundary condition at x = 0, as claimed in [27]. In the resonant case, 
we take an eig enaction,, corresponding to the eigenvalne A To, S efl« := u(l)/u(-l)[ 
and denote by S(9) the free Schrodinger operator restricted to functions in Wf (R \ {0}) 
obeying the interface condition 

/(0+) = 0/(0-), ef(o+) = f(O-); 

then the limit S is given by S(9). 

In fact, resonances in the transmission probability for 5'-like potentials have earlier 
been observed by Christiansen a.o. [12]. In that paper, an exactly solvable model ([ID 
with a specially chosen step function V = a^B was considered. The authors found 
a discrete set of intensities a n (called resonant values) for which partial transmission 
through the limiting <5'-potential occurs; this transmission was shown to rapidly decay 

| Note that in the resonant case any solution of the equation u" = Vu with u'{— 1) = is constant 
outside the convex hull of the support of V; thus the number 9 does not change if in the definition of 
the operator we replace the interval [—1,1] by any interval containing the support of V . 
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as a n becomes larger. The values a n are roots of a transcendent equation depending on 
the regularization of 5'. The findings of [12] are also in contradiction with the results 
of [27] suggesting that the ^'-barrier is completely opaque, no matter what value a 
takes. Exactly solvable models with other piecewise constant potentials as well as non- 
rectangular regular izat ions of 5' have later been studied in [29H3T] . It is also worth 
mentioning the recent results by Zolotaryuk [32H31] . 

In [16] a similar effect was discovered for the family of Schrddinger operators on 
the line of the form 



where ^ G C£°(— 1, 1), W is a real valued potential tending to +oo as \x\ — > oo, and 
a G R is a coupling constant. The map assigning a self-adjoint extension S(a, of the 
operator S' + W to each pair (a, \P) was constructed there. The choice of the extension 
is determined by proximity of the energy levels and the pure states for the Hamiltonians 
with smooth and singular potentials respectively. Two spectral characteristics of the 
profile \l/ are introduced in [TH]: the resonance set E#, which is the a-spectrum of the 
Sturm-Liouville problem —w"+a$fw = on the interval (—1, 1) subject to the boundary 
conditions w'(— 1) = w'(l) = 0, and the coupling function 9^: E^ — > K defined via 
9q,(a) := w a (l)/w a (— 1), where w a is an eigenfunction corresponding to the eigenvalue 
a G E^. In the case when the coupling constant a does not belong to the resonance set, 
S(a, *f>) is just the direct sum of the Schrddinger operators with a potential W on semi- 
axes, subject to the Dirichlet boundary condition at the origin. In the resonant case, 
when a G E#, S(a, acts via S(a, = —f"+Wf on an appropriate set of functions 
obeying the interface condition /(0+) = 9y(a)f(0— ) and 9q,(a)f'(0+) = /'(0— ). 

After we have established the main results of this paper, Prof. Albeverio drew 
our attention to the related work [TJITHEE]- "'" n P^J' ^ e ^ a demonstrated existence of 
'resonant' non-trivial limits for a similar family of the Dirichlet Schrddinger operators 
on the half-line producing in the limit the Robin boundary condition at x = 0. When 
studying the problem of approximating a smooth quantum waveguide with a quantum 
graph, the authors of [UHI] a ls° faced the question on the norm resolvent convergence 
of the family S e of Under the assumptions that V decays exponentially fast at ±oo 
and the mean value j R V of V is non-zero, the authors singled out the set of resonant 
potentials V producing a non-trivial limit of S e in the norm resolvent sense as e — > 0. 
Although the papers [HE] treat a more general situation, the very important case 
j R V = giving S'-type potentials in the limit is excluded from the consideration, and 
the analysis there crucially relies on the fact that j R V ^ 0. In contrast, our approach 
is insensitive to the mean value of V. We also note that our definition of the resonant 
potential agrees with that of the papers [XlfTTj . 

One should keep in mind that the results of the cited papers [HUH EH [16] and of 
this paper do not allow to define the Schrddinger operator with potential 5'. Indeed, the 
limiting operator So is shape-dependent, i.e. it depends on the background potential V, 
so that different approximating families produce different limits. However, these results 




On norm resolvent convergence of Schrddinger operators 



6 



suggest the best choice of an idealized solvable quantum-mechanical model for a realistic 
quantum-mechanical device modelling a dipole. 

We should also note that it is somehow surprising that the limit of S £ exists in the 
norm resolvent sense, although the potentials V £ in general do not converge even in the 
distributional sense. Although V £ {x) — >• for x ^ as e —> 0, the topology of pointwise 
convergence is too weak to force any type of convergence of S e . 

The paper is organized as follows. In the next section, we consider the resonant case 
and prove the convergence result; in section [3], the easier non-resonant case is treated. 
In section HI we specialize the above results to the case of 5'-like potential V = a^> and 
characterize the resonant values of a (i.e., those a for which V is resonant). Finally, 
in the last two sections, we discuss scattering at the resonant values of a and illustrate 
this effect by a simple example. 

Notation. Throughout the paper, VK|(fi), j = 0,1,2, stands for the Sobolev space 
of functions defined on a set Q C M that belong to L 2 (£l) together with their derivatives 
up to order j. The norm in W|(f2) is given by 

ll/llw 2 2 (n) : = (||/"Hi 2 («) + Il/Ill 2 (f7)) 1/2 ) 

where 




is the usual L 2 -norm. We shall write ||/|| instead of ||/||l 2 (m) an d note that H/'Hia^) < 
ll/||w|(n) h y interpolation. 

2. Resonant Case 

In this section, we analyze the more difficult resonant case where A = is an eigenvalue 
of the Sturm-Liouville operator M on (—1, 1) defined in the introduction via ([3]) and 
the Neumann boundary conditions at x = ±1, and denote by u a corresponding 
eigenfunction satisfying the condition u(— 1) = 1. Since the spectrum of Af is simple, 
the function u is uniquely defined. Next, we set 9 := u(l) and denote by S(8) the free 
Schrddinger operator on the line acting via S(6)y = —y" on the domain 

dom S(9) = {ye W*(R\ {0}) | y(0+) = 9y(0-), 9y'(0+) = y'(0-)}. (4) 

The operator S(9) is a symmetric extension by two dimensions of the symmetric 
operator S' with deficiency indices (2, 2) and therefore is self-adjoint. 
Denote by v a solution of the Cauchy problem: 

-v" + V{Z)v = o, £e(-M), v(-l) = o, t/(-l) = i. (5) 

Clearly, we have v G W / 2 2 (— 1, 1); moreover, the Lagrange identity yields (v'u—vu')\_ x = 0, 
i.e., v'(l) = 0- 1 . 
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Lemma 2.1 Fix k 2 E C with Qk 2 ^ 0; then there exists c > with the property that 
for every f E L 2 (IR) and £ > there is q £ E L 2 (M) with \\q £ \\ < ce 1 ! 2 ||/|| such that the 
function y £ = (S(6) — k 2 )~ l f + q £ belongs to domS^ and satisfies the inequality 

\\{S £ -k 2 )y £ -f\\<ce l l 2 H/ll . (6) 

Proof. Fix an arbitrary / E ^(M), set y := {S{9) — k 2 )~ x f and consider the auxiliary 
Cauchy problem 

-< + V(£K = /(e£), ee(-i,i), w £ (-l) = 0, <(-l) = 0. (7) 
We extend i>, w e and the eigenfunction it to the whole line by zero and introduce the 
function 

z £ {x) = (1 - x{x/e))y{x) + y{0-)u{x/e) + ey'{0-)v{x/e)+e 2 w E (x/e), 

with x being the characteristic function of the interval [—1,1]. By construction, the 
function z £ belongs to iy 2 2 (M \ {—e,e}). Although z £ is in general discontinuous at the 
points x = —e and x = e, its jumps and the jumps of its first derivative at these points 
are small. 

To justify this, let [g] a denote the jump of a function g at a point x — a. We first 
observe that (S(0) — k 2 )^ 1 is a bounded operator from L 2 (M) to the domain of S(6) 
equipped with the graph norm; since the latter space is equivalent to W 2 {R\ {0}), there 
exists a constant c\ > independent of / such that 

IMItv|(R\{o}) < ci ll/H ■ (8) 
Now for x = — e we get the estimates 

\[zeU\=W-) - y(s)\< £ \y'(t)\dt < e^Hi/Hw^uo}) < c^*\\f\\, (9) 

HU\=\y'(0-) - y'{-e)\< f \y"{t)\ dt < e 1/2 |blk |(M \ {0} ) < c,e^\\f\\. (10) 

Next, since W|(R \ {0}) C C 1 (1R \ {0}) by the Sobolev embedding theorem, one 
gets a constant c 2 > such that 1 1 ^7 j | c 1 (ir\{o>) < c 2||silw 2 (R\{o}) 101 a U 9 e W 2 (M.\ {0}). 
Using the Sobolev embedding theorem on the interval [—1,1] and properties of solutions 
to the problem (J7J), one arrives at the estimates 

1 1 We 1 1 CI [-1,1] < C2|H|w?(-l,l) < C 3 ||/(£-)||i 2 (-l,l) < C 3 ^ 1/2 ||/|| (11) 

with some constant c 3 independent of /. Therefore, taking into account the equality 
1/(0+) = 9y(0— ), we see that the jump of z £ at x = e can be bounded as 

\[z £ ] £ \=\y(e) - 6y(0-) - ey'(0-)v(l)-e 2 w E (l)\ 

<\y(e)-y(0+)\+e\y'(0-)\\v(l)\+e 2 \w E {l)\ 

<(e 1/2 + c 2 K;(l)|£)||y|| w|(RU0}) +c 3£ 3 / 2 ||/|| <C4e 1/2 ||/|| (12) 

for some C4 independent of / and e E (0, 1). Similarly, keeping in mind that v'(l) = 6~ l 
and y'(0+j = 0~V(O-) and using (EJ) and ([EE]), we get 

|[4] e |=|y'(e) - y'(0-)v'(l)-ew' £ (l)\<\y'(e) - y '(o+)|+e|<(l)| 

< ^llvll**<*\ { o>) W^ll/ll < c 5 s^\\f\\. (13) 
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Figure 1. Functions with the prescribed jumps at the origin 

Let us introduce functions (f and ip as on figure [T]that are smooth outside the origin, 
have compact supports contained in [0, oo), and have the prescribed jumps [<p]o = 1, 
[<//]o = and Mo = 0, fo^o = 1. Set 

GO) = N-e <p{-x - e) - [z' £ ]- E i/;{-x - e) - [z e ] e <f{x - e) - [z' £ } £ ij){x - e); (14) 
then C e = on (-e,e) and, in view of (j9l)- (fTUj) and (fT2]Hfl3lh 

max|Ci fc) (^)| <^ 1/2 ||/|| (15) 

|aj|>e 

for some c 6 > 0, k = 0,1,2, and all e > 0. Clearly, the function y £ := z e + ( £ is 
continuous on R along with its derivative and belongs to W|(R) = domS^. Observe 
that y £ = y for |x| large enough; more exactly, we have 

y e = y + q £ = (S(9)-k 2 )- 1 f + qe (16) 

with 

Q e {x) := -x(x/e)y(x) + j/(0-)u(z/e) + «/(0-)w(x/£)+£ 2 w e (x/£) + Ce(x) 

of compact support. The first four summands above are functions bounded uniformly 
in e and have support contained in [—£,£], and the last one is small due to ( TT5]) . It 
follows from the above estimates that, for a suitable C5 independent of / and e G (0, 1), 

Ikell < cre^ll/H. (17) 
We are now in a position to prove ([6]). If \x\ > e, then 

(S £ - k 2 )y £ = (--^ - k 2 ) (y + q e ) = f- (« - k %. (18) 

If |x| < s, then ( e (x) = 0, and so y £ (x) = z £ (x) = y(0— )«(§) + ey'(0—)v(^)+e 2 w £ (^). 
Thus, 

- k 2 )y £ = 5- 2 y(o-) {-«»(!) + Mf )} + e-V(o-) {-«"(§ ) + v(f Mf )} 

+ {-<(§) + V(§ K(f )} - P X (f )Ve(x) (19) 
= X (f)[/(x)-A; 2 2/£ (x)], 

since both u, v are solutions to equation (jSJ) and w e is a solution to ([7j). Therefore, 

(S e -k 2 )y e = f + r e , (20) 

where 

r £ (x) = -C'J{x) - k 2 [( £ (x) + x(x/e)y E (x/e)} 
= -C(x) - k 2 [q £ {x) + X (x/e)y(x)} . 
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Relations OH]), (EE5J), and ([17]) now yield the estimate 

l|r.|| <ce l/2 \\f\\ (21) 

for a suitable c > c 7 , and (EJ) is proved. □ 



Theorem 2.2 Assume that X = is an eigenvalue of the operator M with an 
eigenfunction u and set 9 := u(l)/u(— 1). T/ien 5 £ converge to S(9) as e in 
the norm resolvent sense. 

Proof. Fix k 2 e C \ K. We conclude from ([20]) and flT£) that 

(5. - = y e - - fc 2 )-v £ = (s(9) - ey'f + (S £ - fc 2 )-V £ 

for each / e L 2 (E). This gives (S e - k 2 )~ x f - (S(9) - k 2 )' 1 f = q £ - (S e - k 2 )~ x r £ , so 
that, by (JI7J and (J2T]), 

IKft " k 2 Y l f - (S(9) - fc 2 )- 1 / II < ||&|| + 11(5. - k 2 y l r £ \\ 

< Hg.ll + I^A; 2 !- 1 !!^!! <c 8 e^ 2 \\f\\ 
for a suitable cs independent of / and s. The theorem is proved. □ 

3. Non-Resonant Case 

Now we study the non-resonant case when A = is not an eigenvalue of N '. Recall 
that SL © S+ denotes the direct sum of the unperturbed half-line Schrddinger operators 
S± = — d 2 /dx 2 on M.± subject to the Dirichlet boundary condition at x = 0; we shall 
prove that S_ © S + is the limit of S £ as e — > 0. 

The above proof still works for the non-resonant case after we have slightly changed 
the corrector q £ of Lemma I2TT1 Let y = (S_ © S + — k 2 )~ l f for some / G L 2 (R) and let 
w be a unique solution of the problem 

-w" + V{£)w = 0, £e(-l,l), w'(-l)=y'(0-), w'(l) = y'(0+). 

It is well known that w obeys the a priori estimate ||w||vi/f < c i(|?/'(0— )| + |y'(0+)|) for 
some constant c\ independent of y'(0±). Combining it with the reasoning of the previous 
section, we find positive constants c 2 , C3, and C4 such that the following inequalities hold: 

IMIc([-i,i]) < c 2 ||w|| W 2 < csllyllwj^o}) < C4||/||. (22) 
We now extend the function w to the whole line by zero and introduce the function 

z £ (x) = (1 - x{x/e))y{x) + ew(x/e)+£ 2 w £ (x/e), 

where w £ is a solution to fl7J) as above. The jumps of z £ and z' £ at the points x = ±e 
converge to as e — > uniformly with respect to the L 2 -norm of /. In fact, taking into 
account the conditions y(0— ) = = y(0+) we obtain 

\[z £ U\ = \ew(-l) - y(-e)\< e\w(-l)\ + J° \y'(t)\dt < c 5 e^ 2 \\fl 

\[z £ ] £ \=\y(e)-ew(l)-e 2 w £ (l)\< f \y\t)\ dt + e\w(l)\+e 2 \w E (l)\ < c^U/U, 

Jo 
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due to (l2"2"j) . Next, the jumps [z' £ ]± e can be estimated as in ( 1TU|) . (ITS]) above. 

Let us introduce the function y £ = z £ + Q, where ( £ is defined as in (JHJ). This 
function belongs to dom S e and can be written in the form 

y £ = {S.®S + -k 2 )~ 1 f + q £ 

with q £ (x) = —x(x/e)y(x) + sw{x / e)-\-e 2 w E {x / e) + ( £ (x). The L 2 -norm of the function 
q £ can be estimated by ce 1 / 2 !!/!!- By calculations similar to those in ( II 8p and (ITS]) , we 
establish that 

(S £ -k 2 )y £ = f + r £ , 

where r £ (x) = —('J(x) — k 2 [q £ (x) + x(x/e)y(x)]. As above, ||r e || < ce 1 / 2 )!/!! for some 
constant c independent of /. 

We now get the following theorem, whose proof is analogous to that of Theorem 12.21 

Theorem 3.1 If X — is not an eigenvalue of the operator Af , then the family S £ of (QJ) 
converges in the norm resolvent sense as e — > to SL © S + . 



4. What is the right Hamiltonian with a ^'-potential? 

Let us consider the formal Schrddinger operators 

d 2 

whose potentials contain the derivative of the Dirac delta function. Here a is a 
strength interaction parameter or a coupling constant taking values in R. Equation 
—v" + a5'(x)v = Xv has no solutions in the space of distributions, except for the trivial 
one. Nevertheless, the formal Hamiltonian H a may be defined in terms of the equation 
— v" = Xv, for i^O, and appropriate boundary conditions at the origin. The main 
question is therefore how to choose these boundary conditions or, in other words, how 
to choose a proper self-adjoint extension of the so-called minimal operator 

S' = ~ t domS' = {ge W 2 (R) | g(0) = g'(0) = 0}. 

The minimal operator S' is symmetric and has deficiency indices (2,2); therefore, its 
self-adjoint extensions form a four-parametric family and there are different possibilities 
to define H a . 

An alternative way is to realize the operator H a as the limit of Hamiltonians with 
regularized potentials 

H e {a^) = ~ + ^l{e- x x) (23) 

with e being a regularization parameter. Suppose that $ G Cq°(R) and supp $ = [—1,1]. 
It is easy to check that e~ 2 ^{e~ l x) — > 5'(x) in the sense of distributions as e — > iff 

/*(0d£ = and / d£ = -1; (24) 
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Figure 2. Odd and more exotic <5'-like potentials 



in this case, we call ^ a 5' -like potential. Plots of some <5'-like potentials are shown on 
figure [2J 

Let us introduce two spectral characteristics of a 5'-\ike potential We denote 
by S^r the set of a-eigenvalues of the problem 



and call E^, the resonant set of It is well known [13] that the a-spectrum of the 
problem ( |25l) is discrete and real and consists of simple non-zero eigenvalues and the 
geometrically simple eigenvalue a = 0; to the latter, there correspond a constant 
eigenfunction and an adjoint function. Also, E# possesses two accumulation points 
at ±oo since the function \1/ changes sign. It is obvious that a belongs to the resonant 
set if the potential a>ty is resonant in the sense of the definition given in the introduction. 
For any non-zero a G E#, we put 



where w a is an eigenfunction corresponding to a. The ratio is properly defined since 
the value w a (—l) is different from 0. Besides, 9q,{a) is always real and does not depend 
on the choice of the eigenfunction. We call 8^: E^ — > R the coupling function of 
Let us introduce the one parameter family of self-adjoint extensions of Sq. 



Hence H(a, is a connected self-adjoint extension with the domain given by (j3J) with 
6 = 8is(a) for the resonant coupling constants a; otherwise, H(a, 1 S r ) is a separated 
extension, namely the direct sum of unperturbed half-line Schrddinger operators subject 
to the Dirichlet boundary conditions. 

Theorem 4.1 Assume \& is a 5 '-like potential. Then the family of Hamiltonians 
H £ (a, converges in the norm resolvent sense as e — > to H(a, 

Proof. This theorem is a reformulation of theorems 12.21 and 13.11 in terms of the poten- 
tial Clearly, A = is an eigenvalue of the operator M with potential V = a^/ if 
and only if the coupling constant a belongs to the resonant set E^. □ 






(25) 
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Note that C^-smoothness of \P is not essential and the theorem also remains valid 
for ^ G Li(R). Moreover, the theorem holds not only for <5'-like potentials but also 
for arbitrary real- valued \P of compact support not necessarily satisfying (I24p . 

Provided that \P is the "shape" of a <5'-like short range potential in the physical 
model under consideration, we can define the limiting Schrodinger operator with a po- 
tential a5' as the operator H(a, \[>). As a result, we find an interesting feature of the 
^'-potentials: the appropriate solvable model is not unique and crucially depends on 
the way in which the derivative of the Dirac delta-function is approximated in the weak 
topology. The shape \1/ is a "hidden" parameter in the conventional formulation of the 
problem on the right definition of the Hamiltonians with <5'-like potentials. 



5. Resonant phenomenon in transmission through a ^'-potential 

It is important to emphasize that the scattering properties of the point dipole also 
depend on the way in which the zero-range limit is realized. We shall show that the 
transmission coefficient depends on the intensity a and the regularizing profile \1/ in such 
a way that for all values of a, the ^'-barrier is completely impenetrable except for the 
set of resonant values, at which there occurs a partial transmission fading away as 
| a | becomes larger. 

First, let us discuss stationary scattering associated with the Hamiltonians H(a, \&) 
and —jgs- We need to consider only the case when the interaction parameter a belongs 
to the resonant set Consider the incoming monochromatic wave e lkx with k > 
coming from the left. Then the corresponding wave function has the form 

Jkx _|_ R e ~ikx if £ < (J, 

Te lkx ifx>0. 



i/j(x, k) 




Here R and T are respectively the reflection and transmission coefficients. As far as 
H(a, ty) = S(0ij,(a)), the matching conditions at the origin 

4>(0+,k) \ = ( 0*{a) \ f if>(0-,k) 
^(0+,k) J \ 6^{aY l ) \ V'(O-fc) 

clearly yield 

' 8*(a) \ / 1 + R 
e^(a)- 1 J I ik(l-R) 

Then one obtains the values of the reflection and transmission coefficients that can be 
expressed via the coupling function as follows (cf. also [1]): 

f 1-6% (a) ( 2Moj 

**{<*) = { l + ltaGL *' T*{a)=\ l + el{a) ltaGL *' (26) 

y —1 otherwise, [ otherwise. 

We emphasize that the reflection and transmission coefficients do not depend on k. 

Next we investigate stationary scattering for the Hamiltonians H e (a, 1 S r ) and — 
and prove that the scattering amplitude converges as e — > to that for the limiting 
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Hamiltonian H(a, The analysis below basically follows a more general approach 
of [23]. 

We are looking for the positive-energy solution of the equation with a ^'-like 
potential 

-y" + ae~ 2 ^{e~ 1 x)y = k 2 y, iel, 
given in the form 

{e lkx + R e~ lkx if x < — e, 

A ek, a) + B v(e~ 1 x, ek, a) if \x\ < e, 

T e ikx if x > e. 

Here u = it(£, x, a) and u = x, a) are the solutions of the equation 

-/ + a*(Ow = A £g(-1,1) (27) 

subject to the initial conditions u{— 1, x, a) = 1, it'(— 1, x, a) = and v (— 1, x, a) = 0, 
v'(—l,x,a) = 1 respectively. The unknown coefficients i?, A, 5, and T can be found 
from the linear system (set x := ek) 



I 



zxe 





1 






1 








it(l, x, a) v{l,K,a) —e 1 



-zxe 



\ / r\ 

A 
B 

\ T ) 



I 



( 



V 



e 

ixe" 





/ 



it'(l,x, a) i/(l, x, a) 

obtained by matching the solution and its first derivative at the points x = ±e. The 
system determinant admits the asymptotic expansion 

A(x,a) = 1/(1,0, a) + ixq(a) + 0(x 2 ), x ->■ 0, (28) 

where q(a) = 2u'(l; 0, a) — it(l; 0, a) — i/(l; 0, a). By Cramer's rule, one obtains 



-u'(l;0,a) + ix[M(l;0,a)-t/(l;0,a)] , 2 n 

-K#(X, a) = r ; p-^r h C(X J, 



it'(l; 0, a) + inq(a) 



T*(x, a) 



— 2ix 



X 



(29) 
(30) 



■u'(l; 0, a) + ixq(a) 

as x -> 0. Here we use the identity it(l, x, a)i/(l, x, a) — u'(l, x, a)i>(l, x, a) = 1 that 
follows from the constancy in £ of the Wronskian of it and v. 

Theorem 5.1 For each k > and a G R i/te reflection and transmission coefficients 
Ry(ek,a) and T^(ek,a) converge towards R^(a) and T^(a) as e — > respectively, 
where the limit values are given by f2t 



Proof. The non-resonant case. Since the equation (1271) for x = coincides with ( 1251 ) 
and a is not a resonant coupling constant, we conclude that it'(l; 0, a) is different 
from 0. From (1291) and (1301) . it immediately follows that Rq,(ek,a) = — 1 + 0(ek) 
and T^(ek, a) = 0(ek) as e — > 0. 

T/ie resonant case. If a G S^, then it'(l; 0, a) = 0. Hence it is an eigenfunction 
of ( 1251) with eigenvalue a. Next, it(l;0, a) = 0#(a) by the definition of the coupling 
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Table 1. Resonant intensities, coupling function and transmission probabilities 



a 


6 9 (a) 


|r*(«)| 2 





1 


1 


18.1747 


-54.9385 


0.00132 


57.1490 


1352.8032 


0.219 ■ 10~ 5 


117.4863 


-32156.4597 


0.387- 10- 8 


199.1756 


755821.4703 


0.704- 10~ n 



function and v'(l;Q,ct) = 9^(a) 1 by the Lagrange identity, which yields the relation 
q(a) = — + 6^(a)). Therefore 



Oil, (a) + 9^(a) 1 + Vy(a) 



T*(ek, a) = - . . j — - + 0(A 2 ) = - 2g *ffi . + 0(e 2 k 2 ) 
K ' 9*(a)- 1 + 6 9 {a) 1 + 91(a) K ' 

as e — > 0, and the proof is complete. □ 



6. An example 

Let us assume that the shape of a short-range potential in an actual model can be 
approximately described as 

{-6x{x + l) if x 6 [-1,0], 
Qx(x-l) ifxe[0,l], 
otherwise. 

The function \& is a 5'- like potential, i.e. it satisfies conditions ( 1241) . As shown above, 
the best choice of the solvable model corresponding to the family of Hamiltonians 

d 2 a T , , . 

in the zero-range limit is given by the operator H(a, Certainly, the resonant set 
and the coupling function 9^ to be found are specific to the given shape 

Since the function \l/ is odd, the set C M is symmetric with respect to the origin, 
namely if a is an eigenvalue of ( 125]) with eigenfunction w a , then —a is also an eigenvalue 
of (125]) with eigenfunction u;_ Q (£) = w a (— £). Furthermore we conclude from this that 

o (- ) 1 

* l ] l«-a(-l) «»a(l) M<*)' 

hence, as follows from ( 126]) . a)| 2 = |i?*(a)| 2 and |T*(— a)\ 2 = \Tq,(a)\ 2 . 

Table [1] lists the first five nonnegative resonant values of a (numerically computed 
using Maple) and the corresponding values of the coupling 9^(a) and the squared 
transmission coefficient |T,j,(a)| 2 . We note that the latter decays very fast and conjecture 
that this will be observed for all £'-like profiles 
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